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Nucleon structure functions can be observed in Deep Inelastic Scattering experiments, but it is 
an outstanding challenge to confront them with fully non-perturbative QCD results. For this pur- 
pose we investigate the product of electromagnetic currents (with large photon momenta) between 
quark states (of low momenta). By means of an Operator Product Expansion the structure func- 
tion can be decomposed into matrix elements of local operators, and Wilson coefficients. For 
consistency both have to be computed non-perturbatively. Here we present precision results for 
a set of Wilson coefficients. They are evaluated from propagators for numerous quark momenta 
on the lattice, where the use of chiral fermions suppresses undesired operator mixing. This over- 
determines the Wilson coefficients, but reliable results can be extracted by means of a Singular 
Value Decomposition. 
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1. Deep Inelastic Scattering 

Historically, Deep Inelastic Scattering gave first evidence that quarks are in fact physical ob- 
jects. More generally it provides insight into the hadron structure functions, and thus into the 
distribution of energy and spin among the hadron constituents (see Ref. |jl|] for a recent review). 

Here we focus on the nucleon structure functions, which can be observed for instance by hard 
leptonic scattering dominated by one-photon exchange, as sketched below. This type of scattering 
only involves one quark, hence chirality is conserved. 




Nucleon Nucleon 



Despite the high energy in this process, perturbative QCD does not lead to fully satisfactory results. 
In particular, power-like IR divergences occur, which imply the notorious renormalon ambiguities 
A consistent evaluation of the nucleon structure function from first principles of QCD has to 
be fully non-perturbative. Hence it is a challenge for lattice studies, and the goal of our project. 

2. Nucleon structure functions on the lattice 

We start from the general ansatz for moments of a lattice nucleon structure function, 

^£(q 2 ) = c (2) (aq)A 2 (a) + < —^-AAa) + ... {higher twists} , (2.1) 

<r 

where a is the lattice spacing, q is the photon momentum being transferred, cW are Wilson coeffi- 
cients and A n are matrix elements (their Lorentz structure is factored out). Traditionally the latter 
are computed by lattice simulations, whereas the coefficients cW are evaluated perturbatively, since 
they only depend on the large photon momentum. However, consistency strictly requires the non- 
perturbative evaluation of both factors [J3j] . Here we present precise numerical results which are 
relevant for c^ . Further steps in this project, leading to ^{q 1 ), are reported in Refs. [Q, ^J]. 

3. Lattice technicalities 

We used configurations that were generated quenched with the Luscher-Weisz gauge action on 
a lattice of size 24 3 x 48 at j3 = 8.45. Based on ro = 0.5 fm, the lattice spacing has been determined 
as a = 0.095(1) fm [0]. We fixed the lattice Landau gauge, which is essential for obtaining finite 
values for the matrix elements; moreover gauge fixing is helpful to reduce the statistical noise. 

Our analysis involves two flavours of degenerate valence quarks, which are represented by 
overlap fermions (we apply an overlap Dirac operator constructed with a Wilson kernel and a 

'The error will be ignored in the following. An alternative method based on F n yielded a ~ 0.105 fm nm. 
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negative mass shift of p = 1.4). The use of chiral quarks suppresses 0(a) lattice artifacts, as well 
as undesired operator mixing. The latter virtue is particularly important in this project; for instance 
it is essential to disentangle the contributions of higher twist, cf. eq. (|Q|). 
Our study includes two bare masses for the degenerate quark flavours, 

m q = 0.01 1 ~ 29 MeV -> m n « 280 MeV and m q = 0.028 ~ 73 MeV -> m n « 440 MeV . 

4. OPE on the lattice 

We use standard lattice electromagnetic currents (it would be computationally expensive to 
work with currents, which are conserved in the framework of overlap quarks). Thus one renormal- 
isation constant (Zy) will be needed The Operator Product Expansion (OPE) decomposes 

the product of two such currents between quark states as follows, 

OPE 

W, v ~ ( ¥ (p)\Mq)jUq)\¥(p))^ I C { ;Jj Mn (q) {wip^L^ip)) (4-D 

m,i,n 

C( m ) : Wilson coefficients, independent of the target and therefore of the quark momentum p 
ff( m > : local operators, relevant to describe the nucleon structure 
\ij : momentum component p^. appears in the operator 

i = 1 . . . 16 : Clifford index , m : index for operators with the same symmetries. 

A truncation of the OPE in eq. ( J4-. lh , and small lattice artifacts, require the scale separation 

p 2 <:q 2 <. (n/af . (4.2) 

Assuming this separation to hold, we consider quark bilinears up to 0(|/?| 3 ), 

y/T>, xfrFD^Y, ijfTD^D^yf , yTD^D^D^y . (4.3) 

The symbol Y captures the full Clifford structure, hence this set includes a frightening number of 
16 • £^ =0 4 d = 1360 operators. However, we choose specific photon momenta of the diagonal form 
q oc (1,1,1,1), which implies a high level of symmetry. To be explicit, we consider three photon 
momenta, 

= \ , I - 2 - 2 GeV i a $ ] = | , k (2) I - 3 - 3 GeV > mf ] = | , k (3) I - 4.4 GeV . 

For q^ and q^ we implement standard boundary conditions (b.c), but qh' is applied along with 
twisted b.c. for the quark fields: in addition to the Euclidean time direction, also two of the spatial 
directions are antiperiodic. This gives access to smaller /^-momenta, which are needed in view of 
condition ((42]), since l^ 1 ) [ is not that large. 

Thanks to our diagonal choice of q, the set of operators reduces to only 67 equivalence classes 
[B]. We classify the corresponding Wilson coefficients according to the number of derivatives in 
the operators that they multiply: 

C\ : no derivative, multiplies C2 ■ ■ Q : one derivative, Bjorken scaling °^\jq 2 

Cj . . .C16 : two derivatives Cn . ..C^ : three derivatives, Bjorken scaling °= 1/ (q 2 ) 2 - 

The coefficients of operators with an even number of derivatives vanish at m q = due to chiral 

symmetry. 
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In each case, our evaluation of C\ . . involves numerous quark momenta p\ , . . . ,pm, see 
Table 0. 
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m q = 0.011 


15 


14 


31 


10 


31 


10 


m q = 0.028 


15 


12 


32 
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31 
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Table 1: The number M of quark momenta, and SV opt of Singular Values (see below) that we used for the 
determination of the Wilson coefficients in each case, i.e. for each quark mass and photon momentum. 

For q( l \ i.e. with twisted b.c, there are less /^-momenta with small p 2 , hence less p-momenta 
are needed for converging results. Thus we measure W^ v — given in eq. (4. 1 ) — off-shell for M = 
15 ... 32 quark momentum sources to determine the Wilson coefficients C\ . . .C(,i- Schematically 
they are given as (the elements and £?jf'' are 4 x 4 matrices capturing the spin components) 



/ W (pi) \ 



V 



\c 67 J 



(4.4) 



Since 16M S> 67 in all our cases, the system is strongly over-determined. Hence we apply 
a Singular Value Decomposition: it selects the n < 67 conditions with "maximal impact" on the 
solution C\ ...C(,j. We order the corresponding Singular Values (SV, analogues to eigenvalues) 
hierarchically. If their magnitude drops rapidly one has favourable conditions to extract a reliable 
result. Fig. [j] illustrates that this is in fact the behaviour that we observed. 



100 



0.01 



m„=0.011 , q <1> , W,,-fit 
q W^-fit 
combined fit 





0.01 



Figure 1: Examples for the rapidly dropping magnitudes of the S V. Due to reflection and rotation symmetries 
some Wilson coefficients coincide theoretically, e.g. the coefficients of ^,-,33 in W33, and of if?, 44 in W44. We 
determine the corresponding SV separately, and by a combined fit, which implements this identity. Left: 
separate and combined SV for m q — 0.01 1, Right: combined SV for both quark masses and q( 2 \ q^\ 

As our next criterion, Fig. || shows how the squared residues in eq. ( fi~4| ) decrease as the number 
of SV involved rises from n = 1 ... 67. 

5. Results for the Wilson coefficients 



Solving eq. (4.4) for C\ . ..Cei employs the inverse SV, so including all of them is not optimal: 
the tiny SV, with large relative noise, tend to distort the result. Therefore we computed the Wilson 
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10 20 30 40 50 60 70 q 2 
n 

Figure 2: Squared residues R 2 for the operator xjffDY (as an example) against the number n of conditions 
(corresponding to the SV). The plots above show results of separate and combined fits (as in Fig. 1), for 
q^ and m q = 0.01 1 (left), and m q — 0.028 (right). The plot below on the left shows the normalised R 2 for 
combined fits at m = 0.028. The final plateau value (at n = 67) is compared in the plot below on the right: 
its decrease for increasing q 2 is roughly linear, so that the relative error remains approximately constant. 

coefficients with a gradually increasing number of SV, n = 1 . . .67; an example is shown in Fig. || 
(left). It displays the most important coefficients, i.e. those of operators with one derivative. The 
only common plateau occurs in the range of 7 ... 13 SV included. To check if this plateau holds 
for all 67 coefficients, we compare the full set obtained with 7, 10 and 13 SV in Fig. ^] (right). We 
observe a striking confirmation of this plateau. The results look similar for other m q and q. The 
optimal number of SV in each case, considering also the impact on <M []|], is displayed in Table ||. 

Next we verify if our number of quark momenta is sufficient: we check if the results change 
significantly as we omit part of them. Fig. ^| shows (with examples) that this is not the case here: 
convergence for an increasing number of /^-momenta is well confirmed. 

In Fig. H (plots above) we compare our results for the Wilson coefficients at to the cor- 
responding values at tree level. As we mentioned in Section 4, C\ (which multiplies 1/^111/0 and 
C-] . . .Ci6 (attached to operators 6^ with y^y v ) vanish at tree level in the chiral limit. We mea- 
sured consistently small values for these coefficients, which indicates that approximate chirality 
and operator mixing are indeed under control, in contrast to previous studies with Wilson fermions 
[B]. Generally the measured Wilson coefficients follow the same pattern as their counterparts on 
tree level, though with significantly reduced absolute values. 

Fig. H (plots below) show the commonly applied Wilson coefficients ratios C m /C^j ee level . The 
dependence on the quark mass is weak, which approves again a safe approximate chirality. On the 
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number of SV 



Figure 3: Results for the Wilson coefficients at m q = 0.01 1 and q^ as a varying number of SV is included. 
Left: the coefficients to the 1 -derivative operators for n = 1 . . .67 reveal a single plateau for n = 7 ... 13. 
Right: the full set C\ ... C(n agrees very well for n = 7, 10 and 13, confirming this plateau generally. 





19 p-momenta 
26 p-momenta 
32 p-momenta 



m g =0.028 , q 
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60 



Figure 4: The Wilson coefficients determined from a varying number of the quark momenta. We show 
results for m q = 0.028, at q^> with 11, 13 or 15 p-momenta (left) and at q^ with 19, 26 or 32 p-momenta 
(right). This hardly changes the results for the Wilson coefficients, hence their convergence appears safe. 

other hand, we observe a strong dependence on the photon momentum q, as expected. Ref. []|] 
discusses the detailed comparison with the theoretically expected Bjorken scaling behaviour (cf. 
Section 4). 



6. Conclusions and outlook 

We have evaluated a set of Wilson coefficients non-perturbatively (partial result were antici- 
pated in Refs. [^]). They refer to twist 2; for the photon momenta that we used, contributions by 
higher twists are suppressed [Sj. The application of chiral lattice quarks has been crucial to control 
the operator mixing. We demonstrated in detail that our results are reliable regarding the number 
of SV and quark momenta included in the evaluation. The measured Wilson coefficients follow the 
pattern of their counterparts at tree level, though with smaller absolute values. 

The structure function ^# (in eq. ([Q|)) is now obtained by means of Nachtmann integration 
over W^v (cf- eq. (4.1)). This is worked out for a single quark in Ref. The final step to a 
fully non-perturbative moment of the nucleon structure function — given by products between the 
matrix elements [[K]] and the Wilson coefficients presented here — is carried out in Refs. [Q, 
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Figure 5: Above: measured Wilson coefficients compared to their tree level values for q( 2 \ Coefficients for 
operators with an even number of derivatives take consistently small values. Below: the ratios C m /C* ee level 
(for non-vanishing denominators). They depend only mildly on the quark mass. 
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